Negative magnetoresistance and phase slip process in superconducting nanowires 
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We argue that the negative magnetoresistance of superconducting nanowires, which was observed 
in recent experiments, can be explained by the influence of the external magnetic field on the critical 
current of the phase slip process. We show that the suppression of the order parameter in the bulk 
superconductors made by an external magnetic field can lead to an enhancement of both the first 
Id and the second I C 2 critical currents of the phase slip process in nanowires. Another mechanism 
of an enhancement of Id can come from decreasing the decay length of the charge imbalance Aq 
at weak magnetic fields because Id is inversely proportional to Aq. The enhancement of the first 
critical current leads to a larger intrinsic dissipation of the phase slip process. It suppresses the rate 
of both the thermo-activated and/or quantum fluctuated phase slips and results in decreasing the 
fluctuated resistance. 

PACS numbers: 74.25.Op, 74.20.De, 73.23.-b 



I. INTRODUCTION 

Recently several experimental groups have observed a 
negative magnetoresistance (NMR) of the superconduct- 
ing nanowires at the temperature lower than the critical 
temperature [J H 0, H B, H ■ In Ref. [IS] the authors 
have directly demonstrated that in their case the effect is 
connected with the suppression of the superconductivity 
in the bulk superconductors caused by the applied mag- 
netic field (it was called as 'anti-proximity effect'). They 
have also found an enhancement of the critical current of 
the nanowires when the magnetic field turns bulk super- 
conductors to the normal state. This result convinced us 
that the observed NMR could be connected with the in- 
crease in the critical currents of the phase slip process. To 
illustrate how the change of the critical current value can 
influence the fluctuated resistance we use a well-known 
model of the point-like Josephson junction with a finite 
capacitance. 

The current-voltage (IV) characteristics of such a junc- 
tion is hysteretic and the parameter which governs the IV 
characteristics is the ratio between the effective "mass" 
(which is proportional to the capacitance C of the junc- 
tion) and the parameter describing the effect of intrinsic 
dissipation (which is inversely proportional to the resis- 
tance R of the junction) Q- In the theory of Joseph- 
son junctions this ratio is called the damping parameter 
f3 c = 2eI c R 2 C/h (7c is the critical current of the Joseph- 
son junction with zero capacitance). The larger [3 C is 
the larger is the hysteresis and the smaller current I r at 
which the voltage vanishes in the junction (see Fig. 1). 
For small values of (3 C the current I r is practically equal 
to I c and the hysteresis is absent. The effect of fluctua- 
tions leads to the appearance of the finite resistance at 
I < I r and the absence of the hysteresis (phase slip pro- 
cess does not stop at I r < I < I c once it is launched by a 



fluctuation) [7( ■ If we increase the value of the current I r 
(by decreasing the resistance of the junction) the fluctu- 
ated resistance at fixed I <C I r decreases (compare gray 
and black dashed curves in Fig. 1). It is the consequence 
of the general rule, that the increase in the intrinsic dis- 
sipation W in the system (in case of Josephson junction 
W = V 2 /R) suppresses both quantum Q and thermo- 
activated fluctuations Q. 
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FIG. 1: Current-voltage characteristics of Josephson junction 
with a finite capacitance and I r \ < I r 2 (parameter j3 c = 5 and 
3 respectively) in the case of zero fluctuations (solid curves). 
Dashed curves show schematically a non-zero voltage response 
at / < I r due to strong fluctuations. In the latter case the 
hysteresis in current- voltage characteristics disappears. 

The aim of the present paper is to show that the mag- 
netic field can enhance the critical currents of the phase 
slip process and hence intensifies the intrinsic dissipation. 
As in the case of Josephson junction it leads to the sup- 
pression of the rate of fluctuations and the decrease of 
the fluctuated resistance. 

However, there is a difference between the results in 
Refs.fl i, 1 i and Refs.|, In Refs. 0, 0, 1 i 
the resistance of superconducting nanowires monotoni- 
cally decreases with the increase of H (up to some critical 
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value) while in Refs. [H, Q the resistance is almost con- 
stant at weak magnetic fields and drops suddenly near 
the critical field of the bulk superconductors. Therefore, 
we conclude that the results are connected with different 
mechanisms of the enhancement of the critical currents 
in the PS process. 

Before considering those mechanisms let us discuss the 
physical meaning of the first and the second critical cur- 
rents of the phase slip process in superconducting wires. 

The paper is organized as follows. In section II we 
observe the critical currents of the phase slip process in 
quasi-lD superconductors. In section III we study the 
influence of the magnetic field on those critical currents. 
Finally, in section IV we discuss our results and make a 
comparison with the experiments and other theoretical 
works. 



II. CRITICAL CURRENTS OF THE PHASE 
SLIP PROCESS 

It has been known for a long time that t he p hase slip 
process in quasi ID superconducting wir es 1 101] is a hys- 
teretic process (see for example review |11| and books 
0, EH ) . If we start from the superconducting state and 
gradually increase the applied current the superconduct- 
ing state becomes unstable at the current l c i (it is ana- 
log of the critical current of Josephson junction I c ) which 
is equal to the product of the depairing current density 
jdep (for wire without defects) on the square of the cross- 
section of the sample S . The periodic in time oscilla- 
tions of the order parameter in one or several points along 
the superconductor destroy the zero-resistance state and 
bring a finite resistance (less than normal one) to the 
system. This process is called a phase slip (PS) pro- 
cess and the points are called as phase slip centers (PSC) 
0, EL E2] • If w e decrease the current below , the phase 
slip process can vanish at the current l c \ < I c % (which 
roughly corresponds to the current I r of Josephson junc- 
tion) . 

The physical origin of the current I c \ was clarified in 
Ref. [13} using the extended quasi-lD time-dependent 
Ginzburg-Landau equations 0, EH ■ It has been found 
that there are two characteristic times which govern the 
dynamic of the order parameter in the phase slip center. 
These are the time relaxation of the absolute value of the 
order parameter t\m and the relaxation time of the phase 
gradient ry^ (which is proportional to the momentum of 
superconducting electrons). 

First characteristic time can be estimated from the 
time-dependent Ginzburg-Landau equation for dynamics 
oflAI 
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|A|(1-|A| 2 -(W) 2 ), 
(1) 



where A = |A|e Ic ^ is the order parameter in Ginzburg- 
Landau equations normalized to its equilibrium value at a 
specific temperature |A| e? = tit B T c v}l 2 (\ - T /Tc) 1 / 2 /n, 



t~gl = ^/(fcs(T c — T)u) is the Ginzburg-Landau re- 
laxation time, £ = (8k B (T c - T)/i:hDy 1 / 2 is a co- 
herence length (D is a diffusion coefficient), 7(T) = 
2\A\ eq (T)TE/fi is the parameter in time-dependent GL 
equations, te is the energy relaxation time for electrons 
near Fermi level and u ~ 5.79 is a number [3, EH- Nu- 
merical analysis shows [l3| that the amplitude of oscil- 
lations of |A| in phase slip center is decreasing with the 
increase of 7 and it is normally much smaller than | A\ eq . 
It allows to neglect the nonlinear term |A| 3 in the right 
hand side (RHS) of the Eq. (1) near the core of the phase 
slip center and it immediately gives us tmi — u~/tgl- We 
can also identify tiai as a relaxation time of the longi- 
tudinal mode in the superconductors tiai — T£fcsT c /A 

urn- 

The second characteristic time Ty^ for long wires L 3> 
Xq (Aq(7 ^> 1) ~ 3> £ is the decay length of the 

charge imbalance [7J) could be estimated by using the 
Ginzburg-Landau equation for dynamics of the phase of 
the order parameter (l3| 



, 2 d 2 ip 



H dej) 



dx 2 
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For Xq > ^ or 7 > 1 the order parameter mainly 
oscillates in a small region (which decreases while in- 
creasing 7) around the phase slip center with the size 
smaller than £. Therefore, we need to estimate ry^, 
in this area. As a result we have ry^, ~ tqlIq^/IXq 
(I is an applied current, Iq — Sh/2eTQL^Pn is propor- 
tional to the depairing current (J C 2 = y / (4/27)/o) J Pn is 
a normal state resistivity and we take into account that 
-d<p/dx{x = 0) = I n (x = 0) Pn /S ~ Ip n /S). 

When T|a| ^$ Ty0 the phase slip process is impossible 
as a periodic one in the time oscillating process [HI at 
/ < I C 2- It allows us to estimate the first critical current 
of long L ^> Lq wires 
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(3) 



where Rps = 2Xqp u /S may be called as a resistance of 
the phase slip process [TtJ • Note, that in case of Joseph- 
son junction current I r is also inversely proportional to 
the intrinsic resistance 0]. 

Due to the above threshold condition there is a voltage 
jump Ay ~ l/ r |A| at I = I c i < I c2 [13]. If current I c \ 
defined by the above expression becomes larger than I C 2 
(at Xq < £ or/and t]a| ^5 t gl) then the voltage gradually 
increases from zero at I = I C 2 and AV = 0. In this limit 
our estimations for tiai and Ty^ become invalid. 

From Eq. (2) it follows that the superconducting elec- 
trons with momentum p ~ V</> being accelerated by 
the gradient of the electrochemical potential fi e ~ch — 
eip + [i c h where p c h might be called as a chemical po- 
tential of the superconducting electrons. In Ref. [16( it 
is shown that —fi c h = eip — p e -ch is also proportional to 
the charge imbalance Q between hole-like and electron- 
like branches of the quasiparticle spectrum in supercon- 
ductors. We may average Eq. (2) over the period of 
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oscillations of | A | and in the case of two-dimensional ge- 
ometry we have 



\ 2 Q AQ~Q = Q. 



(4) 



When the width (w) of the superconducting wire (which 
is connected to bulk superconducting reservoirs) is much 
less than Aq we can leave only the term with the Lapla- 
cian in Eq. (4) near the ends of the wire and solve 
2D Laplace equation. Besides we can neglect the vari- 
ation of Q over the width of the nanowire and solve 
ID variant of Eq. (4) in the wire. As a result we ob- 
tain the charge imbalance at the ends of the nanowire 
Q ~ wQ c /(Xq sui1i(L/2Aq))((5 c is the charge imbalance 
in the phase slip center). Usually w/L <C 1 [1, Q and 
Qo <C Q c even for short wires Aq ^ L. Therefore instead 
of 2D Eq. (4) we may use (in the wire) ID equation 



A; 



d 2 Q 
' dx 2 



-Q = 0, 



(•5) 



with boundary conditions Q(±L/2) = 0, Q(±0) = ±Q C . 

Using Eq. (5) with above boundary conditions it can 
be found that the current I c \ depends on the length of 
the nanowire. If the wire is much longer than the co- 
herence length we may expect that the order param- 
eter distribution in the core of PSC is not influenced 
by the bulk superconductors. Then the dynamics of 
the |A| stays the same as for an infinite wire and both 
T| A | and AV do not suffer any change. From the so- 
lution of Eq. (5) it might be easily seen that the nor- 
mal current in the phase slip center (which is propor- 
tional to the applied one) grows with the decrease of L 
as I n {0) ~ ~dQjdx{x — 0) ~ 1/ tanh(L/2Ag) to provide 
the same charge imbalance Q c ~ AV/2 ~ 1/ T |A| near 
the PSC. Therefore, the shorter the wire is the larger 
is I c i (I c i ~ 1/(Aq tanh(L/2Ag))) and for sufficiently 
short wires it becomes equal to I C 2 . We should note that 
I C 2 does not vary while the wire is much longer than £. 
Therefore, we expect that for short wires L <C Aq the 
hysteresis in current voltage characteristics disappears 
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III. EFFECT OF THE MAGNETIC FIELD 
A. First mechanism 

Let us now discuss how an external magnetic field may 
influence I c \ and I C 2- First mechanism comes from the 
suppression of the order parameter in the bulk supercon- 
ductors. In Fig. 2 we draw the qualitative distribution of 
|A| and Q at H = 0, H < H h c ulk and H > H h c ulk . When 
drawing the curves we have assumed that the NS bound- 
ary forms far from the ends of the wire at H < H c ndk 
and approaches the wire at H > H b c ulk . 

Due to conversing the normal current into the su- 
perconducting one at the NS boundary an additional 
charge imbalance appears at the ends of the wire and 



an effective boundary condition for Eq. (5) becomes 
Q(±L/2) = ±Qo- It brings us the following expression 
for the first critical current 



Ia{H,L) 



-Qo/Qc 



sh(L/2A q ) 



I cl (H = 0,L = oo) 



sinh(L/2AQ) 



(6) 



Current I c i increases with the growth of H because Qo 
changes from zero to the maximal value (with the sign op- 
posite to Q c - see Fig. 2) when the NS boundary touches 
the end of the nanowire (in the latter case the expres- 
sion for I c i was found in Ref. [Hj|). With the increase of 
a magnetic field the point where Q = approaches the 
center of the wire (see Fig. 2). Hence we can say that 
the appearance of NS boundaries effectively shortens the 
superconductor (in sense that 'space' for phase slip pro- 
cess decreases) and it is the reason for an enhancement 
of I cl . 

Because the normal current exists on the finite distance 
from the NS boundaries inside the superconductor the 
current I C 2 is enhanced too (for wires L » Indeed, 
when the normal current penetrates far into the sample it 
decreases the superconducting component of the current 
because I s + I„ = I. Hence, we need a larger applied 
current / to satisfy the condition I C 2 = I s . 
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FIG. 2: Schematic distribution of the charge imbalance and 
order parameter in the superconducting bridge (L — 4Aq) 
with one phase slip center at different values of the magnetic 
field. 

The characteristic length of the discussed mechanism is 
the decay length of the charge imbalance. It means that 
the effect exists only in relatively short wires L < Aq. 
Besides the nanowire should not be wide, otherwise the 
critical field of the wire H c ~ 1/ (£,w) and a bulk super- 
conductor i/^" ifc ~ l/(£ 2 ) become close to each other and 
the magnetic field strongly suppresses |A| in the wire (in 
framework of GL model |A| = (I- [H / H c ) 2 f/ 2 ). It leads 
to increasing t) A | ~ 1/(1 - (H/H c ) 2 and Aq - 1/\A\^ 2 
[IH and hence to decreasing I c i if the effect of the NS 
boundaries is weak. In derivation of Eq. (6) we suppose 
that both T|a| and Aq do not depend on the magnetic 
field. It is true if H <§; H c and at some additional condi- 
tions (see the subsection below). 
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Second mechanism 



The second mechanism of a variation of I c i comes from 
the dependence of pair-breaking mechanisms on the mag- 
netic field due to an orbital effect [ll|[lH ■ In Ref. the 
decreasing Aq with the increase of the applied magnetic 
field was predicted for weak magnetic fields. Because the 
first critical current depends on Aq as I c \ ~ 1/ Aq we can 
expect that I c \ increases in weak magnetic fields. Note 
that in contrast to the first mechanism I C 2 decreases in 
this case. 

The quantitative expression for \q(H) was found in 
Ref. (3 



Aq(#) = 



4Dk B TT E 

n\A\(H) Vl + 7(0)(ff/# c (0)) 2 ' j(T¥ 



(7) 

for temperatures close to T c . The physical reason for the 
dependence of Aq on H is the following. The decay of 
the charge imbalance in superconductors occurs due to 
Andreev reflection process (for quasiparticles with energy 
less than |A| near the NS boundary) or/and due to an 
inelastic electron-phonon interaction (for quasiparticles 
with the energy larger than |A| and along the whole su- 
perconductor). A weak magnetic field almost does not 
influence the order parameter, but it smears the density 
of states of the quasiparticles and makes possible 
Andreev reflection process of the quasiparticles with the 
energy larger than |A|. ft provides a faster relaxation 
of the charge imbalance and decreases Aq. A high mag- 
netic field strongly suppresses the order parameter and 
makes the contribution to Andreev reflections smaller, ft 
increases effective Aq. 

We study this effect in the framework of the time- 
dependent GL equations. We use the field and 
temperature dependent parameter jjj — 7/(1 + 
7(#/ff c ) 2 (l-T/T c ) 1 /2)i/2 inEq. (1,2) and add the term 
{~\A\{H/H C ) 2 ) in the right hand side of Eq . (1) (we 
use dependence H C (T) = H c (0)y/l - T/T c ). For exam- 
ple for tin (7 ~ 100(1 - T/Tc) 1 / 2 ) the effect is expected 
to be weak and for aluminium (7 ~ 10 4 (1 — T/Tc) 1 / 2 ) it 
should be noticeable already at the temperature close to 
T c . 

In our calculations we take T = 0.9T C and several val- 
ues of the parameter 7 (see Fig. 3(a)). At a weak mag- 
netic field I c i increases if the parameter 7 is sufficiently 
high. It occurs not only due to decreasing Aq ~ v/7h 
but also due to decreasing ti^i ~ Jh- The large mag- 
netic field suppresses the order parameter considerably 
and it leads to decreasing l c \ and I C 2 (see the end of 
Sec. Ill A). Because for short wires l c \ does not depend 
on Aq the enhancement of the I c \ due to the discussed 
mechanism should be weakened in such a samples. The 
numerical calculations confirmed that statement (see Fig. 
3(b)). 
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FIG. 3: (a)Dependence of J c i on the magnetic field which was 
found from the numerical solution of Eqs. (1,2) with a field 
dependent parameter j H = 7/(1+7(1?/ H c ) 2 (1-T/T C ) 1/2 ) 1/2 
at different values of 7. (b) Dependence of 7 c i on the magnetic 
field at the fixed value of 7 = 200 and different lengths of the 
superconducting wire. 



IV. DISCUSSION 

We can present simple interpretations of our results. In 
long (L 3> Lq) wires the resistance of the phase slip pro- 
cess Rps is proportional to 2Aqp„/5 fl71 |. In short (L < 
Lq) wires Rps is proportional to the length of the wire 
Rps — Lp n /S because the electric field is close to zero 
in bulk superconductors. When a short wire is bounded 
by the normal metal the resistance of the phase slip pro- 
cess decreases by the value which is proportional to the 
resistance of the NS boundaries Rps — Lp n /S — 2Rns- 
When one applies a magnetic field Rps can decrease due 
to decreasing Aq (in long wires) or due to the appearance 
of Rns ( m short wires) and as in the case of Josephson 
junction current I c \ ~ AV/Rps increases. The intrinsic 
dissipation W ~ V 2 /Rps increases too and it suppresses 
the fluctuations of the phase of the order parameter and 
results in decreasing the fluctuated resistance. 

The extended time-dependent Ginzburg-Landau equa- 
tions are valid for a narrow temperature interval close 
to the critical temperature \T C — T\ < h/ksTE where 
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£(T) > Le = \JDte- However, we may suppose that 
traces of the high temperature dynamics should exist 
at low temperatures too. It was experimentally found 
for a tin nanowire with L = 6fj,m the presence of sev- 
eral phase slip centers (see Fig. 3b in Ref. [2(|) at 
T = Q.5K ~ 0.12T C . First phase slip brings a finite resis- 
tance which is equivalent to the resistance of the piece of 
the wire with the length about 1.5 \xm and hence using 
SBT theory [l7| we obtain Xq(Sti) ~ 750nm. It is much 
larger than the coherence length in tin at this tempera- 
ture (£(Sn,T = 0) < 55nm). 

In Ref. [2lJ the S-behavior of the current voltage 
characteristics of superconducting nanowires was found 
in a voltage driven regime at low temperatures. Time- 
dependent Ginzburg-Landau equations with a large value 
of 7 give qualitatively the same result [H, [2l| . The volt- 
age jump AV was extracted from the experimental data 
and showed a qualitative agreement with the theoretical 
temperature dependence AV(T) (l3j . 

In the recent paper [22j the pronounced hysteresis of 
the IV characteristic of the Pb nanowire with several 
phase slip centers was observed at a low temperature. 
In the preceding experiments [H, [2l| the hysteresis was 
hidden by a strong external noise. Proper filtering sup- 
pressed the noise [22| and revealed the hysteretic behav- 
ior with well identified currents I c \ and I C 2- 

The above experiments support the idea that the low- 
temperature properties of the phase slip process resemble 
such ones at high temperatures. Therefore, we expect 
that our results based on the numerical solution of Eqs. 
(1,2) and Eqs. (4,5) are applicable for qualitative analysis 
of the phase slip process at low temperatures. 

The magnetic field dependence of the first and second 
mechanisms of the current enhancement are rather differ- 
ent. In the case of the first mechanism the enhancement 
occurs at the moment when bulk superconductors switch 
to the normal state. In the second mechanism the current 
Id increases smoothly with H (see Fig. 3). Therefore, 
we believe that the first mechanism may be responsible 
for the effects observed in Refs.Jl, Q while the second 
mechanism is connected with the experiments in Refs. 

Indeed, the theory for the first mechanism clarifies 
why the 'anti-proximity effect' 0, Q is weakened in long 
Zn nanowires with L = SOfim (Aq(Zti) ~ 22/im) and 
is absent in Sn nanowires with lengths L — 6 — 30/im 
(Xq(Sti) ~ 750nm - see our estimation above). They 
seem to be too long to observe the 'anti- proximity effect'. 

In Ref. [|| direct measurements demonstrated an 
increase of the critical current in the superconducting 
nanowires. We identify this current as the first critical 
current of the phase slip process. The second critical cur- 
rent was not observed in the experiment due to a high 
rate of fluctuated PSC which reveals itself in the finite 
resistance of the wire at I — ► 0. The absence of hysteresis 
is also typical for Josephson junction with a high rate of 
fluctuations Q. 

To observe the first mechanism of NMR (due to cur- 



rent enhancement) the width w or diameter d of the 
wire should be relatively small to provide a condition 
H c ~ l/(w, d) > H h c ulk . Otherwise the magnetic field 
H = H^ ulk strongly suppresses the order parameter in 
the wire. It leads to decreasing both critical currents J c i, 
I C 2 (see and of Sec. Ill A) and results in a negligible effect 
of the normal boundaries. We believe that it is the rea- 
son for the observed dependence of the 'anti-proximity 
effect' on the width of the nanowires in Ref. [3]. The 
authors found no effect for a wire with d = 70nm ~ £/2 
and pronounced effect for a nanowire with d — 40 ~ £/4. 




(H/H) 2 

FIG. 4: Dependence of the order parameter averaged over the 
width of the film on the applied magnetic field for different 
values of w (the result is obtained in framework of Ginzburg- 
Landau equations) . It shows a relatively stronger suppression 
of I A in films with w > 2£ than in narrow films w < 2£ for 
the same ratio H/H c . 

The maximal enhancement of the current I c \ does not 
depend on w, d (while w, d < £) in the second mech- 
anism. The invariance is connected with a universal 
dependence of the order parameter on a magnetic field 
A(H) = A(0)(1 - (H/Hc) 2 ) 1/2 for narrow samples. For 
wider samples w > £ the order parameter decreases faster 
with the growing ratio of H/H c (see Fig. 4) and the en- 
hancement of I c i becomes weaker. For very wide samples 
the vortices enter the sample at H <C H c ~ H C 2 and that 
suppresses the order parameter even stronger. The con- 
clusion is that the second mechanism of enhancement of 
I c \ is maximal in narrow quasi-ID samples with w^d < £. 

We can speculate that the observed in Ref. [l[ sup- 
pression of the NMR for wide bridges is connected with 
exceeding width of the sample over the coherence length. 
For example the minimal NMR occurs for a sample with 
w = lOOnm (see Fig. 4 of Ref. |l|). Taking into ac- 
count the coherence length of bulk lead (£(0) = 83nm) 
and a dirty limit expression for the coherence length 
£ - y/hD/\A\ we have £(T = 1.6K) ~ 50nm (for sam- 
ples with T c ~ 2.5K). Therefore, for the widest sample 
our rough estimation gives w ~ 2£. Besides from our 
theory of the second mechanism of NMR it follows that 
the maximal current enhancement (or maximal suppres- 
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sion of the resistance) occurs at H* ~ H c ~ l/(w,d) for 
w ,d < (see Fig. 3a). In Refs. [2, H[ qualitatively the 
same dependencies were observed. 

Our both mechanisms give a suppression of NMR at 
approaching T c (similar to experiments [H, 0, H, B @] ) • 
The first mechanism becomes noticeable if the length of 
the wire is not very short. Otherwise the current Z c i is 
equal to J C 2 and I C 2 decreases for wires with L < 4£(T) 
due to a strong proximity effect from the normal banks. 

The second mechanism of NMR is effective only when 
the term 7(0)(i? / 'H c (0)) 2 is larger than the unity (see 
Eq. (7)). It is obvious that for some temperatures close 
to T c the above condition fails because H < H C (T) <C 
H c (0). Besides with increasing temperature the charge 
imbalance length increases as Aq ~ Aq(0)(1 — T/T c ) -1 / 2 
and the length of the sample becomes smaller than 
Aq. It also suppresses the second mechanism of NMR 
(see Fig. 3b). 

The proposed mechanisms of the negative magnetore- 
sistance are different from those studied in Refs. [23ll24j. 
The authors of work [23| supposed that additional resis- 
tance due to NS boundaries stabilizes the superconduct- 
ing phase. In our approach the intrinsic dissipation grows 
due to decreased intrinsic resistance of the phase slip pro- 
cess and it suppresses the fluctuations in the system. 

In Ref. [24| a new channel of dissipation in supercon- 
ducting wires was proposed which can be suppressed by 
an external magnetic field. It would be interesting to 
compare the contributions of that channel and the sec- 
ond mechanism of NMR studied in the present paper. 

Our mechanisms of current enhancement are also 



rather different from a critical current enhancement pre- 
dicted in Refs. [26|, [2?], HI]. In those works the en- 
hancement of current I C 2 (using our terminology) was 
found due to suppression by applied magnetic field the 
pair-brea king resulting from total 'spin-flip'-|- 'non-spin- 
flip ' rate [28|]. This process leads to decreasing the cur- 
rent I c \ because Aq grows according to that effect. In 
our second mechanism the increased pair-breaking (due 
to orbital effect) decreases current I c i but enhances the 
current l c \. 

Our first mechanism has rather different behavior on 
magnetic field and cannot be confused with theory of 
Refs. [13, HI] . The direct way to distinguish among the 
second mechanism of a current enhancement is predicted 
in our paper and Refs. (27l . HH is to study the samples 
of different lengths and widths (with other close param- 
eters). In contrast to Refs. [27l l28j the current enhance- 
ment in our second mechanism depends on the length of a 
nanowire and becomes weaker for wires with the length 
L < 2Aq. The second difference is that the maximal 
current enhancement does not depend on the width of a 
wire for narrow samples w < £ while in (2rl I27] ] the strong 
dependence on w was predicted even for such a narrow 
wires. 
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